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EVALUATION  OF  SOLAR  FRARES  AND  ELECTRON  PRECIPITATION 
BY  NITRATE  DISTRIBUTION  IN  ANTARCTICA 


Summary 

During  the  beginning  of  the  funding  period,  most  of  the  time 
devoted  to  project  research  was  spent  in  Antarctica  (November 
1990  to  mid  January  1991).  A  firn  core  was  drilled  by  hand  to  a 
depth  of  29  meters  at  Windless  Bight  on  the  Ross  Ice  Shelf. 

This  dill  site  was  located  about  10  km  from  that  of  the  1988-89 
core.  The  main  result  of  the  study  is  that  all  of  the  major 
peaks  identified  as  resulting  from  ionization  caused  by  SPEs  in 
1972,  1946,  and  1928  that  were  found  in  the  1988-89  core  could 
also  be  identified  in  the  analytical  sequence  from  the  1990-91 
core.  Dating  of  the  new  core  is  based  mainly  on  the  snowfall 
record  from  nearby  McMurdo  Station,  snow  stratigraphy  and  high 
resolution  conductivity  data.  Following  the  Antarctic  field 
season,  we  were  able  to  obtain  a  set  of  snow  samples  that  had 
been  collected  by  the  Chinese  member  of  the  International  Trans- 
Antarctica  Expedition  that  crossed  the  entire  continent  on  foot 
in  a  7  month  period  in  1989-90.  The  analysis  of  these  95 
samples  showed  nitrate  flux  that  correlates  closely  with  known 
spatial  distribution  of  electron  precipitation  in  the  south 
polar  region.  A  new  apparatus  has  been  built  for  field  analysis 
on  a  continuous  basis  of  nitrate  and  conductivity  in  a  melt 
stream  derived  from  the  vertical  melting  of  ice  cores.  Such  a 
system  will  provide  an  improved  method  for  dating  the  individual 
ice  layers,  as  well  as  maximum  accuracy  and  precise  matching 
between  the  nitrate  and  conductivity  profiles. 
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EVALUATION  OF  SOLAR  FRARES  AND  ELECTRON  PRECIPITATION 
BY  NITRATE  DISTRIBUTION  IN  ANTARCTICA 


I.  Distribution  of  Nitrate  Content  in  the  Surface  Snow  of 
the  Antarctic  Ice  Sheet  Along  the  Route  of  the  1990 
International  Trans-Antarctica  Expedition 


Previous  work  showed  that  nitrate  measured  at  very  high 
resolution  (1.5  CM)  in  snow  depositional  sequences  in  Antarctica 
could  be  correlated  with  short  term  phenomena  such  as  solar 
proton  events  (Dreschhoff  and  Zeller,  1990).  It  was  clear  that 
deposition  of  the  ionization  products  in  the  snow  is  strongly 
dependent  upon  precipitation  and  atmospheric  conditions  during 
and  immediately  after  the  event.  Information  about  the 
geographic  distribution  of  the  nitrate  fallout  over  Antarctica 
was  limited  to  only  a  few  sites  however.  A  unique  opportunity 
to  examine  this  aspect  of  the  nitrate  distribution  and  to  test 
the  hypothesis  that  atmospheric  ionization  from  solar  charged 
particles  is  responsible  for  a  significant  portion  of  nitrate 
was  presented  to  us  by  a  set  of  surface  snow  samples  collected 
by  the  International  Trans-Antarctica  Expedition  foot  traverse. 
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Figure  1.  Nitrate  tlux  in  surface  sanples  across  the  Antarctic 
continent  sanplod  during  the  International  Trans-Antarctica 
Expedition. 
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The  set  of  95  samples  of  the  upper  25  cm  was  collected  by  Qin 
Dahe  of  the  Chinese  Academy  of  Science  at  roughly  equal 
distances  along  the  5736  kilometer  route  from  27  July,  1989  to 
3  March,  1990.  Samples  are  distributed  along  a  track  from 
65”05"S,  59'’35'W,  through  90*  S,  to  66'’33'S,  95°39'E,  which 
represents  geomagnetic  latitudes  ~50  S,  West  longitude,  to 
~77*S,  East  longitude.  The  profiles  of  nitrate  concentration 
and  flux  along  the  route  were  plotted  and  indicate  that 
(especially  at  the  higher  elevation  of  the  polar  plateau)  the 
distribution  is,  at  least  in  part,  controlled  by  the 
precipitation  of  electrons.  The  highest  flux  was  found  in  areas 
where  electron  precipitation  is  known  to  be  unusually  high 
(Sheldon,  et  al.,  1987).  These  results  are  shown  in  Figure  1, 
and  details  of  this  work  are  presented  in  Appendix  A.  Further 
interpretation  of  this  data  is  anticipated.  As  pointed  out  by 
D.  Smart  and  M.  A.  Shea  (personal  communication,  1991),  data 
such  as  the  auroral  boundary  index  compiled  by  Madden  and 
Gussenhoven  (1990)  could  result  in  a  more  comprehensive 
explanation  of  variation  in  that  portion  of  the  nitrate  profile 
that  traverses  the  area  surrounding  the  South  Pole  Station.  A 
preliminary  discussion  of  the  nitrate  distribution  relative  to 
electron  precipitation  in  the  south  polar  region  will  be 
presented  at  the  AGU  Fall  meeting,  San  Francisco,  1991  (See 
Appendix  B) . 

II.  Results  from  Windless  Bight,  Antarctica  on  Major  Solar 
Flares  during  Solar  Cycles  22  to  14. 

During  the  1990-91  field  season  in  Antarctica  we  drilled  29 
meters  of  firn  core  by  hand  and  made  simultaneous  analyses  for 
both  nitrate  and  liquid  conductivity  using  a  Beckman  Model  160 
UV  spectrophotometer  and  an  Orion  Model  160  conductivity  meter. 
The  detector  cells  were  arranged  in  sequence  so  that  flow  from 
the  UV  absorption  cell  was  discharged  directly  intc  the 
conductivity  cell  and  readouts  were  obtained  simultaneously. 

The  seqence  for  nitrate  could  be  compared  with  that  obtained 
during  the  1988/89  Antarctic  field  season  in  which  a  62  year 
high  resolution  analytical  sequence  from  the  Ross  Ice  Shelf 
showed  strong  anaomalies  that  could  be  directly  associated  with 
specific  manjor  solar  flare  events  (Dreschhoff  and  Zeller, 

1990).  This  profile.  Figure  2a,  represented  cycles  16  through 
21.  Large  solar  proton  events  (SPEs)  were  easily  detectable  in 
cores  sampled  continuously  at  high  resolution  (1.5  cm)  and  these 
primary  peaks  were  over  6  standard  deviations  above  the  mean  (in 
Figure  2  the  peaks  are  labeled  A,  B  and  C) .  At  this  resolution, 
secondary  peaks  caused  by  meteorological  effects  associated  with 
individual  snow  storms  (year  1988)  also  occurred  at  intensity 
levels  up  to  a  maximum  of  4  standard  deviations  above  the  mean. 
It  was  found  that  identifiable  peaks  in  the  nitrate 
concentration  are  critically  dependant  on  factors  such  as  the 
timing  of  the  event  relative  to  the  polar  winter  night  and  the 
local  meteorology  that  influences  snow  accumulation. 
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Figure  3.  Nitrate  concentrations  in  relative  units  (UV 
absorbance)  at  two  drill  sites  at  Mindless  Bight,  for  Antarctic 
season  1988-89  (2a)  and  Antarctic  field  season  1990-91  (2b).  A 
new  Becksan  Model  160  which  is  less  sensitive  by  about  30  I  was 
used  in  1990-91.  The  two  sequences  are  separated  by  -  lo  )«  and 
represent  -  21  seters  of  firn  core  end  a  tine  period  of  -  62 
years.  Local  seteorology  and  ice  flow  produce  differences  in 
snow  deposition  and  snow  cospaction  changes  progressively 
through  tine  resulting  in  differences  in  thickness  of  yearly 
layers.  Nevertheless,  the  prinary  concentration  peaks 
associated  with  flares  in  1972,  1946,  and  1938  are  clearly 
visible  (labeled  as  A,  B  and  C).  Dating  has  been  accospllshed 
by  snowfall  records,  visible  snow  stratigraphy,  D/H  variations 
and  conductivity. 


Furthermore,  after  deposition,  diagenetic  processes  related  to 
the  compaction  of  the  snow  tend  to  distort  the  shape  of  the 
original  nitrate  signal  by  diffusion  within  ice  crystals  and  by 
mixing  that  occurs  as  the  crystals  grow  together  in  the 
compaction  process. 

These  processes  introduce  a  randomization  factor  that 
increases  through  time  and  plastic  flow  of  the  ice.  Nitrates 
scavenged  from  the  atmosphere  (mostly  troposphere)  are  thought 
to  reside  on  the  outer  surfaces  of  the  snow  crystals  (W.  H. 
Zoller,  Personal  communication,  1991).  If  this  is  true,  it 


might  be  expected  that  nitrate  from  this  source  would  be 
subject  to  major  redistribution  during  the  recrystallization 
processes  that  occurs  in  surface  and  near  surface  snow 
associated  with  compaction.  On  the  other  hand,  nitrate  from 
stratospheric  sources  is  probably  incorporated  into  snow 
crystals  by  processes  that  cause  a  major  difference  in 
distribution  within  the  snow  crystals. 

Polar  stratospheric  clouds  (PSCs),  in  some  cases,  consist 
of  very  high  concentrations  of  nitric  acid  (McElroy  et  al., 
1988).  We  hypothesize  that  precipitating  snow  crystals  may 
contain  nucleation  centers  derived  from  PSC  particles.  Once 
they  enter  the  troposphere,  these  nucleation  centers  are 
surrounded  by  condensing  ice  and  may  survive  the  diagenetic 
process  without  significant  redistribution.  For  this  reason, 
nitrate  peaks  originating  from  PSCs  or  other  sources  in  the 
stratosphere,  like  those  identified  as  A,  B,  and  C,  are  more 
likely  to  retain  their  identity.  Some  preliminary  evidence  for 
this  hypothesis  is  presented  in  Figure  2b.  This  is  the  data 
series  acquired  during  1990/91  Antarctic  field  season  equivalent 
in  depth  to  the  record  from  1988-89.  The  complete  record  is 
shown  in  Figure  3a.  It  is  important  to  recognize  that  these 


Fiqure  3.  Conplete  29  nater  firn  core  froa  tha  1990-91  field 
season  showinij  nitrato  conountrutions  in  relative  units  of  UV 
absorbance  (3a)  and  liquid  conductivity  (3b).  Headouts  of 
nitrate  and  conductivity  were  obtained  siBultaneously  on  the 
sane  sanplc.  Arrow  indicates  position  of  dust  layer. 
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primary  peaks  are  present  in  both  sequences  in  spite  of  the  10 
kilometers  distance  that  separates  the  drill  sites  and  the 
differences  in  meteorological  conditions  and  ice  flow  patterns 
in  the  ice  shelf.  Nevertheless,  differences  in  the  expression 
of  the  secondary  peaks  are  clearly  visible. 

During  the  field  collection  and  analysis  of  the  1990-91 
firn  core,  we  found  a  visible  pale  brownish  layer  approximately 
3  cm  thick,  near  the  bottom  of  the  core  at  a  depth  of  27.8 
meters.  The  layer  was  especially  notable  in  the  diffused 
lighting  present  in  the  trench  at  the  drill  site.  This  layer 
was  analyzed  in  the  normal  sequence  of  core  processing  and  was 
found  to  contain  extremely  high  levels  of  nitrate  and  also  to 
show  conductivities  that  were  very  high  (Figure  3b.).  Using  a 
petrographic  microscope,  we  found  that  the  ice  contained 
particulates  that  included  highly  angular  fragments  of  glassy 
material  as  well  as  some  crystalline  minerals  (see  Figure  4.). 

In  this  case,  it  is  most  likely  that  the  high  nitrate  and  high 
conductivity  of  the  samples  is  related  to  the  deposition  of  a 
layer  of  volcanic  dust  on  the  snow  surface.  Some  of  the  nitrate 
contained  in  this  layer  could  have  been  derived  from  chemical 
processes  taking  place  in  the  eruption  cloud  and  from 
electrostatic  discharges  which  are  known  to  accompany  volcanic 
eruptions.  Furthermore,  if  the  dust  were  deposited  on  the  snow 
surface  during  summer,  it  would  have  the  effect  of  increasing 
absorption  of  solar  radiation  and  would  have  greatly  increased 


I 

I 

I 
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Figure  4.  Photomicrograph  of  dust  contained  in  ice  samples  from 
the  dust  layer  shown  in  Figure  3.  Haqnif ication  is  X-200  in 
cross  polarized  light  showing  angular  glassy  aggregates  and 
birefringent  crystalline  grains.  Most  probably  derived  from  the 
volcano. 
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evaporation,  causing  concentration  of  soluble  constituents  in 
the  surface  layer.  In  any  sequence  of  data  obtained  from  ice 
cores,  it  is  necessary  to  be  able  to  distinguish  between  large 
primary  nitrate  peaks  that  are  the  result  of  giant  solar  flares 
and  other  large  peaks  that  result  from  different  phenomena  like 
volcanic  eruptions  or  major  dust  storms.  In  this  case  the 
sequence  from  the  core  collected  in  1990-91  at  Windless  Bight 
provides  an  illustration  of  this  type  of  problem  (see  Figure 
3a.).  Details  of  the  1991  data  sequence  will  be  discussed  at 
the  AGU  Fall  Meeting  in  San  Francisco,  1991  (see  Appendix  C). 
This  portion  of  the  research  was  funded  in  part  by  grant 
NSF/DPP  8919190. 

111.  Apparatus  Development  for  the  1992  Greenland  Field  Season 

Recently  new  apparatus  has  been  designed  and  built  by  us  and 
the  Center  for  Reasearch  Inc.,  University  of  Kansas  which  will 
allow  ultra-high  resolution  (-2  mm)  measurement  of  nitrates  and 
conductivity  based  on  an  automated  continuous  analytical 
system.  This  apparatus  called  the  Polar  Ice  Core  Analysis 
System  (PICAS),  has  been  constructed  and  is  now  under  test  in 
our  laboratory.  The  equipment  permits  a  melt  stream  from  a 
heated  probe  that  will  melt  vertically  through  the  ice  core  to 
be  pumped  by  a  micropump  through  the  well-tested  Beckman  Model 
160  uv  spectrophotometer  analyzer  cell.  The  effluent  from  this 
cell  will  then  be  passed  in  turn  through  the  microcell  of  the 
Orion  Model  160  conductivity  apparatus  as  was  done  in  the  system 
that  we  used  in  Antarctica  in  1990-91.  Readout  from  both 
instruments  is  fed  to  a  computer  in  the  control  unit  where  the 
data  can  be  displayed  for  preliminary  evaluation  and  then  stored 
on  disk.  Construction  of  the  melter  probe  has  not  been 
completed  but  a  preliminary  design  has  been  tested.  The 
micropump  has  been  acquired  and  mounted  on  a  variable  speed 
electric  motor  and  is  now  undergoing  testing.  The  pump  and  all 
major  components  are  stainless  steel.  The  integrated  unit  will 
be  field  tested  in  Colorado  this  winter.  The  entire  system 
weighs  less  than  100  lbs  and  is  designed  for  field  operation  in 
polar  climates. 

After  the  field  test  and  any  redesign  that  may  be 
necessary,  this  system  will  be  used  during  the  1992  Greenland 
field  season.  Our  team  (Dreschhoff  and  Zeller)  is  scheduled  to 
work  at  the  Summit  site  (see  map  in  Figure  5.)  which  is  the  deep 
drilling  site  (GISP  II)  for  the  United  States  and  European 
glaciology  research  projects.  The  Polar  Ice  Coring  Office 
(PICO)  will  drill  a  100-120  meter  deep  core  for  us  (see  letter. 
Appendix  D.)  We  expect  to  make  all  analyses  of  samples  at  the 
drill  site.  It  would  be  impossible  to  complete  analysis  in  the 
field  in  a  single  season  without  the  PICAS  system. 
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Figure  i.  Map  ot  crccnlanU  showing  the  location  of  the  tield 
area  at  the  iiuaait  drill  site.  . 
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V.  Interactions  and  Activities 

1 .  )  Based  on  our  experience  in  using  a  carbon  dioxide  laser  for 
cutting  ice  cores,  we  have  also  designed  a  laser  ice  core  drill 
which  is  described  in  detail  in  Appendix  E. 

2. )  November  1990  -  January  18,  1991.  Antarctic  field  season 
and  acquisition  >  80  years  of  high  resolution  nitrate  and 
conductivity  data.  We  also  comtacted  Prof.  Westphal  in 
Antarctica  where  we  discussed  neutrino  detection  systems  based 
on  the  use  of  polar  ice  as  a  Cerenkov  medium.  We  had 
previously  suggested  this  to  Prof  Halzen  of  the  University  of 
Wisconsin  in  1988.  (see  acknowledgements  in  Appendix  F.) 

3. )  May  1991  was  spent  at  Zhengzhou  University  Physics 
Department  as  part  of  the  World  Bank  -  Provincial  Universities 
Development  Project. 

4. )  From  June  1991  to  September  1991  our  laboratory  was  visited 
by  Prof.  Cin  Oahe,  Lanzhou  Institute  of  Geology  and  Geocryology 
Chinese  Academy  of  Sciences,  China. 

5. )  We  have  contacted  Dr.  Julie  Palais  of  NSF  Division  of  Polar 
Programs  to  exchange  data  on  the  dust  layer  found  near  the 
battom  of  our  firn  core. 

6. )  Dr.  Joan  Feynman  visited  the  Physics  Department  at  the 
University  of  Kansas  in  August  1991  and  subsequently  spent  one 
day  visiting  our  lab.  A  statistical  comparison  between  our 
nitrate  profiles  was  discussed. 
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Apendix  A 

Preprint  | 
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Abstract 

Previous  work  showed  that  nitrate  measured  at  very  high  resolution  (1,5  cm)  in  snow 
depositional  sequences  in  Antarctica  could  be  corrected  with  short  term  phenomena 
such  as  solar  proton  events  (Dreschhoff  and  Zeller,  1990).  It  was  clear  that  deposition  of 
the  ionization  products  in  the  snow  is  strongly  dependent  upon  precipitation  and 
atmospheric  conditions  during  and  immediately  after  the  event.  Information  about  the 
geographic  distribution  of  the  nitrate  fallout  over  Antarctica  was  limited  to  only  a  few 
sites  however.  A  unique  opportunity  to  examine  this  aspect  of  the  nitrate  distribution 
and  to  test  more  fully  the  hypothesis  that  atmospheric  ionization  from  solar  charged 
particles  is  responsible  for  a  significant  portion  of  nitrate  was  presented  to  us  by  a  set  of 
surface  snow  samples  collected  by  the  International  Trans-Antarctica  Expedition  foot 
traverse.  The  set  of  95  samples  of  the  upper  25  cm  was  collected  by  one  of  us,  (Qin),  at 
roughly  equal  distances  along  the  5736  kilometer  route  from  27th  July,  1989  to  3rd 
March,  1990.  Samples  are  distributed  along  a  track  from  65°05'S,  59®35'W,  through  90°S, 
to  66°33'S,  95'’39'E,  which  represents  geomagnetic  latitudes  50°S,West  longitude,  to  77°S, 
East  longitude.  The  profiles  of  nitrate  concentration  and  flux  along  the  route  were 
plotted  and  indicate  that  (  especially  at  the  higher  elevation  of  the  polar  plateau)  the 
distribution  is,  at  least  in  part,  controlled  by  the  interaction  between  the  magnetosphere 
and  the  flux  occurs  in  areas  where  electron  precipitation  is  known  to  be  unusually  high 
(Sheldon  et  al.,  1988) 


Introduction 
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The  nitrate  content  in  the  snow  and  ice  at  several  sites  in  Antarctica  have  been 
reported  by  previous  authors  (Paker  et  al.,  1978;  Zeller  and  Parker,  1981;  Delmas, 
1986;  Legrand,  1987;  Dreschhoff  and  Zeller,  1990;  Mayewski  and  Legrand,  1990). 
These  previous  studies  focused  not  only  on  the  nitrate  concentration,  and  the 
distribution  and  temporal  variations,  but  also  on  the  sources  about  which  the 
authors  have  different  opinions.  However,  these  glaciological  investigations, 
including  the  nitrates,  in  the  Antarctic  ice  sheet  are  limited  to  a  few  sites,  such  as 
South  Pole,  Vostok  Station,  Dome  C,  Ross  Ice  Shelf,  etc.  No  systematic  investigation 
of  the  nitrate  content  has  been  made  in  the  surface  snow  of  the  major  geographic 
zones  of  the  Antarctic  ice  sheet,  especially  in  a  single  austral  field  season,  although ' 
major  differences  in  concentration  at  sites  were  reported  as  being  related  to  their 
location  with  reference  to  the  auroral  zone  (Zeller  an  Parker,  1981).  This  paper 
presents  the  results  of  the  nitrate  analysis  of  surface  snow  samples  collected  along 
the  1990  International  Trans-Antarctica  Expedition  (1990  ITAE)  route  from  about 
50°S  magnetic  latitude.  West  longitude,  through  90®  S  to  about  77°  S  magnetic 
latitude.  East  longitude. 


Sampling  and  analytical  procedures 

The  1990  ITAE  completed  the  crossing  of  the  Antarctic  continent  along  the  longest 
route,  starting  from  Seal  Nunataks  (65°05'S,  59°35'W)  on  27th  July,  1989,  and 
arriving  at  Mirnyy  Station  (66°33'S,  95°39'E)  on  3rd  March,  1990.  The  traverse 
covered  the  northern  part  of  the  Larsen  Ice  Shelf  and  the  southern  part  of  the 
Antarctic  Peninsula,  Palmer  Land,  Ellsworth  Mts.,  South  Pole  (90°00'S),  Vostok 
(78°28'S,  106°48'E)  and  Komsomolskaya  Stations  (74°05'S,  97°27'E)  (Fig.  1).  The  total 
distance  of  5736  km  was  completed  within  220  days.  95  samples  of  25  cm  depth  of 
surface  snow  were  collected  in  polyethylene  bottles  along  this  route  at  equidistant 
points.  The  bottles  were  new  and  were  washed  three  times  using  Milli-Q-Water  to 
clean  each  bottle  and  permit  us  to  maintain  accuracy  in  the  ppb  range  for  ionic 
measurement  of  the  snow  samples.  The  washed  bottles  were  packed  in  clean  bags  in 
the  clean  room  and  were  opened  only  during  sampling  for  a  short  time.  Specific 
steps  were  adopted  to  keep  the  samples  clean  during  sampling,  for  example,  to  face 
the  windward  direction,  to  wear  gloves  and  to  clean  the  surface  of  the  profile  of  the 
snow  pit,  so  that  contamination  was  avoided. 
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All  of  the  samples  were  kept  frozen  until  processed.  The  nitrate  content  of  the 
samples  was  analyzed  at  the  Space  Technology  Center,  University  of  Kansas,  USA, 
by  a  Beckman  160  UV  Selectable  Wavelength  Spectrophotometer.  Absorbance  is 
determined  with  an  analytical  precision  of  <  2%,  and  when  calibrated  against 
standards,  indicates  that  the  mean  concentration  for  the  95  snow  samples  collected 
on  the  ITAE  was  22.2  ppb  ±1  ppb  (NO  ‘  -  N). 

Before  the  processing  of  the  samples  in  Kansas,  these  95  samples  had  been 
analyzed  for  12  ions  by  ion  chromatography  in  the  Laboratory  de  Glaciologie  et 
Geophysique  de  L'Environnement,  France.  All  of  the  samples  have  been  melted  in 
the  clean  room  in  France  for  3-6  hours,  then  immediately  refrozen  and  kept  in  cold 
storage  for  shipment  to  the  United  States.  For  the  ionic  analyses,  all  samples  were 
handled  in  a  manner  that  avoided  any  contamination,  so  that  the  accuracy  of  the 
next  analyses  would  not  be  compromised  (Qin  Dahe  et  al.,  in  preparation).  The 
excellent  agreement  between  the  analytical  results  in  France  and  those  at  the 
University  of  Kansas  suggest  strongly  that  these  precautions  were  successful. 

Results 

Table  1  presents  the  results  of  the  nitrate  analysis  of  the  samples,  as  well  as  the 
detailed  basic  data  of  the  sampling.  In  Table  1  the  sampling  date,  the  position  and 
the  distance  from  Station  1  (Seal  Nunataks)  are  recorded.  The  elevation  value  is 
obtained  from  the  large  contour  map  of  Antarctica  edited  by  Drewry  (1982), 
according  to  the  records  of  station  position.  The  mean  annual  accumulation  rate  has 
been  collected  from  previous  publications,  and  its  source  is  indicated  in  Table  1,  and 
the  annual  nitrate  flux  is  the  product  of  the  measured  nitrate  concentration  and  the 
snow  accumulation  rate.  The  data  on  the  accumulation  rate  is  available  only  for  61 
stations,  for  which  most  of  the  data  possess  high  accuracy,  especially  in  the  segment 
from  Station  47  to  85  (Table  1).  Assuming  the  mean  density  of  the  25  cm  surface 
snow  is  0.3  g/ cm3  (Qin  Dahe  and  Ren  Jiawen,  1991;  Lipekove,  1980),  the  25  cm  snow 
would  cover  roughly  1-3  average  year's  snow  falls  in  most  of  the  61  sampling 
stations.  Taking  into  account  the  complex  nature  of  snow  accumulation,  it  can  be 
argued  that  the  25  cm  surface  snow  was  deposited  in  the  recent  past.  That  means  the 
nitrate  flux  in  the  25  cm  snow  represents  present  nitrate  deposition,  in  recent  snow 
accumulation. 
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Fig.  2  and  Fig.  3  show  the  distributions  of  the  nitrate  concentration  at  95  stations 
and  of  nitrate  flux  at  61  stations  in  recent  snow  along  the  5736  km  route  across 
Antarctica. 


Discussion 

The  distribution  of  the  nitrate  concentration  in  25  cm  surface  snow  along  the  route 
of  the  1990  ITAE  in  Fig.  2  shows  that  the  concentration  in  the  inland  segment  is 
higher  than  that  in  the  coastal  segments.  To  be  precise,  the  nitrate  concentration  in 
the  surface  snow  of  the  Antarctic  Plateau,  from  Station  47  to  Station  84,  where  the ' 
elevation  is  between  2600  m  and  3560  m  a.s.l.,  is  generally  higher  than  that  at  the 
other  stations,  where  the  sampling  sites  are  limited  to  the  Larsen  Ice  Shelf,  Antarctic 
Peninsula  or  to  localities  within  a  600  km  distance  from  the  coast  lines.  Fig.  2  also 
shows  that  the  highest  concentration  values  are  not  in  the  center  of  the  Plateau,  but 
in  the  segments  on  the  route  of  2  to  3  degrees  of  latitude  on  both  sides  of  the  South 
Pole.  However,  the  nitrate  concentration  in  the  25  cm  surface  snow  alone  can  not 
describe  quatitatively  the  overall  characteristics  of  nitrate  distribution  in  Antarctica. 
In  Fig  3.,  the  relationship  between  the  nitrate  flux  and  the  route  distance  is  plotted 
and  it  provides  an  obvious  illustration  of  the  nitrate  deposition  along  the  route.  In 
Fig.  3,  there  are  only  61  nitrate  flux  values,  because  it  is  impossible  to  collect  the 
complete  accumulation  data  of  all  stations.  This  figure  shows  that  the  nitrate  flux 
value  is  higher  along  the  coast  or  in  the  near-coast  regions,  and  generally  lower  in 
the  inland  region.  However,  with  the  exception  of  three  stations  located  outside  of 
the  inland  segment,  the  highest  nitrate  flux  value  is  at  Station  48  (88°38'S,  96°26'W). 
In  Fig.  3,  we  also  can  observe  that  the  lower  nitrate  flux  values  are  concentrated  in 
the  center  of  the  Antarctic  Plateau  segment,  mainly  from  Station  60  (85°33'S, 
105‘’40'E)  to  Station  80  (74°44'S,  98°41'E).  In  this  segment  the  nitrate  concentrations 
shown  in  Fig.  2  are  high  and  show  strong  variations  from  station  to  station. 

To  account  for  the  above-mentioned  distribution  we  must  consider  the 
mechanisms  by  which  the  nitrate  is  produced  and  transported  to  Antarctic  surface 
snow.  In  ,  the  prior  studies,  Zeller  and  Parker  (1981)  had  considered  a  total  of  14 
possible  mechanisms  for  the  origin  of  nitrate  in  Antarctic  snow  and  ice.  Clausen  and 
Langway  (1989)  also  indicated  that  the  nitrate  impurity  observed  in  the  Antarctic 
snow  has  a  multiple  source  origin.  From  the  earlier  analyses,  the  nitrate  sources 


13 


f 


could  include  continental,  atmospheric,  stratospheric  contributions,  and 
anthropogenic  activity  effects,  consisting  of  biomass  burning,  soils  and  dust, 
oxidation  of  NH^,  photochemical  processes,  lighting,  and  several  effects  of  cosmic 
related  factors.  But  Prospero  and  Savoie  (1989)  analyzed  the  observed  data  of  nitrate 
collected  from  a  Pacific  island  network,  and  found  the  nitrate  concentration  in  the 
air  of  the  Northern  Pacific  to  be  three  times  greater  than  that  of  the  Southern  Pacific. 
They  found  in  the  Northern  Pacific,  that  the  higher  nitrate  concentration  is 
correlated  with  Asian  dust  transport,  and  by  contrast,  suggested  that  the  nitrate 
concentration  in  the  Southern  Pacific  region  is  related  to  oceanic  'background' 
sources.  Petit  et  al.  (1991)  analyzed  the  excess  deuterium  in  recent  Antarctic  snow' 
and  concluded  that  the  moisture  of  the  precipitation  in  the  center  of  the  Antarctic 
ice  sheet  is  derived  from  the  warmer  ocean  of  the  southern  lower-mid  latitude.  In 
the  coastal  region  of  the  Antarctic  ice  sheet,  the  moisture  of  the  precipitation  is 
derived  from  the  adjacent  oceans.  From  these  observations,  it  follows  that  biomass 
burning,  continental  and  anthropogenic  sources,  can  basically  be  neglected  as 
substantial  contributors  to  nitrate  in  Antarctic  snow  and  ice. 

Because  of  the  remoteness  of  Antarctica  from  other  continents,  upper  atmospheric 
sources  make  an  important  contribution  to  the  nitrate  contained  in  Antarctic  snow 
and  ice.  This  interpretation  furnishes  a  partial  explanation  for  the  nitrate  flux 
distribution.  Our  investigation  shows  that  the  nitrate  flux  in  the  inland  segment  is 
lower  than  that  in  the  coastal  or  near-coastal  segments  of  Antarctica,  especially  at  the 
edge  of  the  ice  sheet  (Fig.  3).  Accumulation  rates  are  very  low  in  the  interior  of 
Antarctica,  and  the  moisture,  as  well  as  some  of  the  nitrate,  derives  from  the  lower 
latitude  Southern  Ocean  where  the  'background"  of  nitrate  is  actually  quite  low 
(Prospero  and  Savoie,  1989).  But  in  the  coastal  region,  the  snow  accumulation  is 
very  high,  with  greater  nitrate  flux  from  the  local  air  masses.  Even  so,  the  nitrate 
concentration  remains  quite  low  (Table  1  and  Fig.  2).  In  the  coastal  region,  there  is  a 
large  extent  of  exposed  soils  and  rocks  in  moraines,  nunataks  and  mountains. 
Claridge  and  Campbell  (1968)  showed  that  soil  in  the  ice-free  areas  is  frequently  very 
high  in  nitrate.  It  appears  likely  that  the  atmospheric  contribution  determines  the 
basic  nitrate  flux  distribution  in  the  surface  snow  of  the  Antarctic  ice  sheet  especially 
in  the  inland  region  as  well  as  in  the  coastal  region.  The  area  of  the  immediate  coast 
is  probably  influenced  strongly  by  the  locally  exposed  soils  which  contribute 
additional  nitrate. 
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However,  our  data  have  revealed  some  unexpected  results  shown  in  Fig.  2  and 
Fig.  3.  Specifically,  the  highest  nitrate  concentrations  are  located  on  both  sides  of  2-3 
degrees  of  latitude  from  the  South  Pole  (Fig.  2).  The  highest  nitrate  flux  in  the 
inland  segment  is  also  located  in  this  area  (Fig.  3).  These  samples  were  collected 
during  the  Antarctic  midsummer  period,  1989  (Table  1).  Moreover,  in  Fig.  3  high 
nitrate  flux  also  appeared  around  78°S  latitude  in  the  Ellsworth  Mountain  region. 
These  samples  were  collected  in  late  October,  1989.  It  appears  that  in  these  segments, 
we  have  to  be  concerned  about  three  component  contributions  to  the  spatial 
distribution  of  nitrate  concentration  and  flux  across  the  continent.  First,  a 
background  component  is  subject  to  general  circulation  and  storm  system  transport 
over  the  Antarctic  Ice  Sheet.  Second,  superimposed  anomalies  in  nitrate  flux  can 
result  from  ioni2ation  in  the  auroral  zone  and  thirdly  in  the  lower  ionosphere  (<  D- 
region)  caused  mostly  by  highly  penetrating  electron  precipitation.  Several  authors 
have  discussed  the  changes  in  the  chemical  composition  of  the  upper  and  middle 
atmosphere  due  to  ionization  by  charged  particles  (Swider,  1977;  Jackman  et  al.,  1980; 
Thone,  1980;  Rusch  et  al,,  1981;  Solomon  et  al,,  1981;  Garcia  et  al.,  1984) 

If  the  nitrate  anomalies  can  be  accounted  for  by  atmospheric  ionization  products, 
the  nitrate  fallout  at  any  specific  collection  site  must  be  related,  in  part,  to  the 
location  of  the  site  relative  to  the  geomagnetic  latitude  (see  Fig.  4).  Specifically,  an 
area  about  2  to  3  degrees  latitude  from  the  geographic  South  Pole  (geomagnetic 
latitude  75  degrees  and  Station  51  in  the  ITAE  traverse)  should  be  expected  to  exhibit 
auroral  zone  contributions  when  compared  to  ITAE  stations  between  South  Pole 
and  Vostok.  Lanzerotti  et  al.  (1984)  and  Rosenberg  et  al.  (1987)  pointed  out  that  the 
conjugation  of  the  magnetic  field  lines  and  spatial  extent  of  open  or  closed  magnetic 
field  lines  mostly  above  70  degrees  geomagnetic  latitude  are  a  strong  function  of 
geomagnetic  local  time  and  the  variation  in  flux  of  charged  particles.  We  would 
expect  that  any  contributions  of  nitrate  as  an  upper  atmospheric  ionization  product 
should  be  at  a  minimum  in  the  very  high  geomagnetic  positions  along  the  ITAE 
route.  This  coincides  with  results  from  the  all-sky  images  taken  at  South  Pole 
Station  (Mende  and  Rairden,  1987).  In  addition,  auroral  absorption  or  ionospheric 
absorption  of  radio  waves  by  electrons  that  penetrate  deeply  into  the  lower 
ionosphere  has  been  found  to  exhibit  a  maximum  at  about  66  degrees  magnetic 
latitude  (Krishnaswamy  and  Rosenberg,  1985). 
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These  types  of  results  by  riometer  measurements  coincide  with  rocket 
observations  and  model  calculations  on  the  coupling  of  relativistic  magnetospheric 
electrons  to  the  middle  atmosphere  (Baker  et  al.,  1987).  In  addition,  there  is  a  strong 
longitudinal  variation  in  the  penetrating  electron  precipitation  in  the  Southern 
Hemisphere  (Sheldon  et  al.,  1987).  In  fact,  Sheldon  et  al.  (1988)  estimated  that  in 
certain  areas  of  Antarctica,  the  rate  of  electron  precipitation  and  therefore  rate  of 
ionization  that  could  cause  changes  in  the  chemical  composition,  of  the  atmosphere 
would  be  particularly  high  because  of  its  association  with  the  magnetic  South 
Atlantic  Anomaly.  ITAE  station  27  (Table  1  and  Fig.  4)  would  clearly  fall  into  this 
zone. 

Some  of  our  nitrate  values  from  three  sations  along  the  ITAE  can  be  compared 
with  previous  results  reported  in  the  literature.  Herron  (1982),  Legrand  and  Delmas 
(1986)  and  Legrand  et  al.  (1984)  calculated  the  nitrate  (in  ppb)  fluxes  at  Siple  (some 
hundred  years  in  the  Holocene),  South  Pole  (1959-1969)  and  Vostok  (Holocene)  to  be 

11,  7.6  and  1.3  kg/km2  a(N02'),  respectively.  Our  results  at  the  same  sites  are  1.8, 1.5 

and  0.6  kg/km2.a(NO  ”  -  N),  respectively.  The  difference  between  the  fluxes  at  these 
sites  may  be  due  to  the  fact  the  nitrate  data  reported  here  represent  samples  from  the 
most  recent  time  period  of  about  two  years,  i.e.  1988-1989.  However  the  general 
characteristic  of  these  two  groups  of  results  shows  that  the  highest  value  is  at  Siple, 
the  lowest  at  Vostok,  with  South  Pole  showing  an  intermediate  flux.  This 
distribution  of  nitrate  fluxes  does  not  contradict  the  apparent  relationship  between 
the  collection  sites  and  geomagnetic  latitudes,  including  the  effects  of  the  position  of 
the  auroral  zone.  For  example,  due  to  the  fact  that  Vostok  Station  is  closest  to  90 
degrees  geomagnetic  latitude,  contributions  from  the  auroral  zone  should  be 
comparatively  low  throughout  all  seasons.  (Mende  and  Rairden,  1987). 

A  significant  difference  between  nitrate  data  reported  in  the  literature  and  the 
results  of  our  measurements  is  to  be  found  in  the  flux  ratio  between  Siple  and 
Vostok  Stations.  Data  from  previous  studies  indicates  a  ratio  of  8.5  whereas  our 
measurements  show  only  3.  Without  more  samples  and  data  relating  to  the  exact 
nature  of  the  snow  sequence  sampled  by  the  previous  workers  at  these  collection 
sites,  these  differences  are  difficult  to  evaluate  and  they  may  actually  represent  a 
change  in  nitrate  deposition  that  has  taken  place  in  the  recent  past.  We  know  that 
deposition  from  polar  stratospheric  clouds  containing  nitric  acid  is  possible  (McElroy 
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et  al.,  1988;  Wofsy  et  al.,  1990)  and  that  there  is  a  depletion  of  stratospheric  ozone 
over  Antarctica  (Fiocco  et  al.,  1989;  Mayewski  and  Legrand,  1990).  We  also  know  that 
there  is  an  interannual  variability  in  the  intensity,  position  and  location  of  the  polar 
vortex  over  the  Antarctic  Ice  Sheet  (Bowman,  1986).  In  short,  the  differences  in  the 
nitrate  flux  ratios  between  Siple  and  Vostok  may  be  the  result  of  significant  changes 
in  ozone  depletion  as  well  as  solar  variability  over  Antarctica  over  the  past  few 
years. 


Conclusions 

In  this  paper  the  nitrate  concentration  and  the  flux  in  the  surface  snow  of 
Antarctica  along  the  route  of  the  ITAE,  5736  km,  has  been  presented.  The 
distributional  characteristics  of  the  nitrate  has  been  defined  by  our  results  for  this 
profile  accross  the  remote  polar  continent  which  is  separated  from  other  land 
masses  by  the  bulk  of  the  Southern  Ocean.  Three  mechanisms  seem  to  be  important 
in  determining  the  spatial  distribution  of  nitrate  in  the  top  25  cm  of  the  surface 
snow  in  the  regions  traversed  by  the  ITAE.  These  three  mechanisms  are  related  to 
the  Southern  Ocean,  auroral  activity  as  well  as  anomalous  electron  precipitation 
and  the  processes  related  to  ozone  depletion.  Although  short-lived  effects  from  solar 
proton  events  undoubtably  are  present  they  could  only  be  identified  as  anomalies  in 
high-resolution  samples  of  snow  sequences  extending  over  at  least  one  year 
(Dreschhoff  and  Zeller,  1990). 

The  local  coastal  ocean  and  the  remote  southern  lower -mid  latitude  ocean  provide 
a  partial  origin  for  the  nitrate  transported  by  atmospheric  circulations  and  deposited 
in  the  Antarctic  snow  and  ice  in  the  coastal  and  inland  regions.  The  variation  in  the 
precipitation  rate  in  Antarctica  causes  the  basic  distribution  of  nitrate  flux,  the 
higher  flux  in  the  coastal  region  and  the  lower  in  the  inland.  For  this  reason,  the 
Southern  Ocean  furnishes  the  largest  contribution  to  the  nitrate  budget  in  Antarctic 
snow  and  ice.  In  some  segments  of  our  sampling  route  the  nitrate  flux  is  strikingly 
increased  in  zones  of  auroral  activity  and  enhanced  electron  precipitation.  These 
phenomena  do  not  seem  to  be  simple  coincidence.  Upper  atmosphere  physics 
studies  provide  extensive  evidence  to  support  this  interpretation. 

While  we  believe  that  we  have  distinct  evidence  of  a  magnetospheric  effect  on  the 
nitrate  fallout  along  the  ITAE  traverse,  the  investigation  presented  here  is  based 
upon  a  small  number  of  samples  (95)  collected  over  a  short  time  (July  1989  to  March 
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1990)  during  a  period  of  extreme  variability  in  charged  particle  emission  from  the 
sun.  Clearly,  to  substantiate  this  conclusion,  future  investigations  will  have  to  pay 
special  attention  to  magnetospheric  effects  as  well  as  spatial  and  temporal  location  of 
the  collection  sites. 
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Figure  Captions 


Fig.  1.  Antarctic  map  showing  the  route  of  the  1990  International  Trans- Antarctica 
Expedition  (ITAE)  and  the  sampling  stations. 

Fig.  2.  Nitrate  concentrations  in  25  cm  depth  of  surface  snow  on  the  Antarctic  Ice 
Sheet  and  elevation  along  the  route  of  the  ITAE. 

Fig.  3.  Nitrate  flux  in  25  cm  surface  snow  and  elevation  along  the  route  of  the  ITAE. 
Fig.  4.  Map  showing  the  ITAE  traverse  relative  to  the  geomagnetic  latitudes  and 
nominal  average  position  of  the  auroral  zone  (shaded  area)  for  geomagnetic  local 
noon  at  the  South  Pole  (adapted  from  Weller  et  al.,  1987). 
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Table  1 


Table  1.  The  results  of  nitrate  concentration  of  95  samples  and  the  interrelated  data. 
The  references  for  accumulation  rate  are  as  follows:  (1).  Kotlyakov  et  al.,  1977;  (2). 
Giovinetto  and  Bull,  1987;  (3).  Peel  and  Clausen,  1982;  (4).  Hamilton  et  al.,  1971;  (5). 
Shimizu,  1964;  (6).  Taylor,  1965;  (7).  Jouzel  et  al.,  1983;  (8).  Jouzel  et  al.,  1990;  (9). 
Lipenkov,  1980. 
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Apendix  B 


Slectron  Precipitation  and  Nitrate  Deposition 
Along  the  Route  of  the  1990  International 
Trans-Antarctic  Expedition 

Jin  Dahe  (Lanzhou  Inst-  of  Glaciology  and  Geo¬ 
cryology,  Chinese  Academy  of  Sciences, 
Lanzhou,  Gansu  730000,  China) 

2.  J.  Zeller  and  G.  Dreschhoff  (Both  at;  Space 
Technology  Center,  2291  Irving  Hill  Drv. 
Lawrence,  KS  66045) 

A  set  of  surface  snow  samples  collected  by 
:he  International  Trans-Antarctica  foot 
iraverse  constitutes  a  unique  opportunity  to 
:est  the  hypothesis  that  ionization  from  solar 
:harged  particles  is  responsible  for  a 
significant  portion  of  the  nitrate  in  Antarctic 
snow.  Ninety  five  samples  of  the  upper  25  cm 
^ere  collected  by  one  of  us  (Qin)  at  roughly 
jqual  distances  along  the  5736  km  route  from  27 
ruly,  1989  to  3  March,  1990.  These  samples  are 
iistributed  along  a  track  from  65‘’o5'S,  59^35'W 
Longitude,  through  90®S,  to  66^33 'S,  95*^39 'E, 
[/hich  represents  geomagnetic  latitudes  ~50^S,  W 
:o  ~77®S,  E  longitude.  Nitrate  concentrations 
md  flux  along  the  route  are  plotted  relative 
:o  distance  and  the  altitude  curve  is 
superimposed  on  the  plots.  Altitude  clearly 
iffects  precipitation  rates  and  influences 
:he  fallout  of  trace  compounds  to  the  surface. 
:n  addition,  the  results  show  that  the  nitrate 
:lux  across  the  Antarctic  continent  is,  at 
.east  in  part,  controlled  by  the  interaction 
)etween  the  magnetosphere  and  the  upper 
itmosphere.  The  results  also  show  that  the 
lighest  nitrate  flux  occurs  in  an  area 
rhere  electron  precipitation  is  known  to  be 
inusually  high. 


1.  1991  Fall  Meeting 


2.  000901157 


3.  (a)  G.  Dreschhoff 
Space  Tech  Center 
2291  Irving  Hill  Dr. 
Lawrence,  KS  66045 
(b)  913-864-4517 

4.  U 

5.  (a)  U03  Paleoastro- 
physics  and  natural 

variation  of  cosmo- 
genic  isotopes 
(b)  2475  Polar 
ionosphere 
2716  Energetic 
particles,  precip. 

7.  0% 

8.  $50  check  incl. 

9.  C 
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Apendix  C  | 


Ionization  Products  from  Solar  Flares 
Retained  in  Antarctic  Snow 

E.  J.  Zeller  and  G.  Dreschhoff  (Both  at:  Space 
Technology  Center,  2291  Irving  Hill  Drv. 
Lawrence,  KS  66045) 

J.  Feynman  (Jet  Propulsion  Laboratory,  4800 
Oak  Grove  Drv.  Pasadena,  CA  91109) 

Over  a  seven  year  period,  high  resolution 
chemical  analysis  for  nitrate  ion  in  firn 
sequences  from  the  Windless  Bight  area  of  the 
Ross  Ice  Shelf  in  Antarctica  have  been  made. 
These  time  series  have  been  assembled  and 
confirm  that  nitrate  concentrations  in  polar 
ice  and  snow  can  be  used  to  provide  reliable 
signals  of  past  solar  flare  activity.  It  has 
been  found  that  major  solar  proton  events 
produce  strong  anomalies  that  are  tracable  from 
nitrate  sequences  separated  by  as  much  as  10 
kilometers.  On  the  other  hand,  anomalies  that 
are  related  to  individual  snow  storms  and 
depositional  conditions  occur  only  in  a  single 
time  series  and  appear  to  be  determined  by  the 
specific  conditions  at  the  site.  These  data 
present  a  basis  for  quantitiative  evaluation  of 
the  nitrate  anomalies  versus  the  proton  flux 
from  the  associated  SPE's.  A  companion  paper 
(Qin,  Zeller  and  Dreschhoff)  shows  that 
electron  precipitation  associated  with  the 
auroral  zone  may  be  responsible  for  the  spatial 
distribution  of  nitrate  fallout  over  the 
Antarctic  continent  and  plays  a  major  role  in 
determining  the  nitrate  background  over  broad 
areas. 


1.  1991  Fall  Meeting 


2.  000900622 

3.  (a)  E.  J.  Zeller 
Space  Tech  Center 
2291  Irving  Hill  Dr. 
Lawrence,  KS  66045 
(b)  913-864-4516 

4.  U 

5.  (a)  U03  Paleoastro- 
physics  and  natural 
variation  of  cosmo¬ 
gonic  isotopes 

(b)  7514  Energetic 
particles 
7519  Flares  and 
mass  ejections 

7.  25%  Solar  Physics 
127,333,1990 

8.  $50  Check  incl. 
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NATIONAL  SCIENCE  FOUNDATION  Apendix  D1 

1800  G  STREET.  N.W. 

WASHINGTON.  D.C.  20550 

Dr.  Gisela  Dreschhoff 
Radiation  Physics  Lab 
University  of  Kansas  Space 
Technology  Center 
2291  Irving  Hill  Drive 
Lawrence,  Kansas  66045-2969 


Dear  Gisela,  September,  13  1991 

Thanks  for  reminding  me  that  I  still  had  not  responded  to  your 
letter  of  July  8,  1991  regarding  your  request  for  assistance  from 
NSF  for  transport  and  shallow  drilling  at  Summit,  Greenland  during 
the  1992  field  season.  Please  contact  PICO  and  ask  them  to  come 
up  with  a  budget  for  the  support  you  require  (separating  the 
drilling,  camp  support  (if  any)  and  transportation).  As  I 
understand  it,  the  Air  Force  will  pick  up  your  transportation  costs 
directly  by  making  arrangements  with  the  109th.  I  suppose  the 
flights  on  which  you  and  your  cargo  are  transported  to  GISP2  will 
cost  PICO/NSF  less  than  a  similar  flight  with  just  PICO/NSF  cargo 
and  personnel.  Eventually  we  can  discuss  with  the  Division  of 
Grants  and  Contracts  (DGC)  the  exact  mechanism  by  which  you  will 
transfer  funds  from  your  Air  Force  grant  to  DPP/NSF  for  the 
drilling  support  you  get  from  PICO. 

In  the  near  future  someone  from  NSF  will  be  sending  you  a 
questionaire  to  fill  out  so  that  we  may  include  your  program  in  our 
application  to  the  Danish  Commision  for  permission  to  work  in 
Greenland.  Chuck  Myers  in  DPP  usually  handles  those 
questionaires.  If  you  do  not  receive  one  by  late  Fall  let  me  know 
so  that  I  can  get  him  to  send  you  one.  Finally,  as  soon  as  next 
summer's  schedule  becomes  clear  we  can  discuss  when  you  might  go 
up  to  Summit.  I  suspect  the  month  of  August  will  be  the  best  for 
us  but  we  can  discuss  that  at  a  later  time. 

Let  me  know  if  you  have  anymore  questions  about  this  matter. 
Good  luck  with  your  work  and  please  say  hello  to  Ed  for  me. 


Sincerely, 


Polar 


M.  Palais 
Glaciology 


% 


cc;  Dr.  J.  Kelley,  PICO 

Dr.  H.  Zimmerman,  DPP/NSF 
Dr.  Ted  DeLaca,  DPP/NSF 
Dr.  Dennis  Peacock,  DPP/NSF 
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Apendix  El 


Development  of  Laser  Ice  Cutting  Apparatus 

Edward  Zeller,  Gisela  Dreschhoff,  Claude  M.  Laird 

Radiation  Physics  Laboratory, 

Space  Technology  Center 
University  of  Kansas, 

Lawrence,  Kansas,  66045 

During  the  1990/91  field  season  at  Windless  Bight  near 
Ross  Island,  our  team  introduced  the  use  of  a  25  watt 
continuous  infrared  carbon  dioxide  laser  as  a  field  device  to 
cut  individual  firn  cores  for  sample  preparation*  The  test  was 
successful  and  permited  this  device  to  be  employed  on  a  routine 
basis  in  field  operations,  carbon  dioxide  laser  cutting 
systems  have  the  advantage  that  the  beam  is  emitted  at  an 
infrared  wavelength  that  is  absorbed  in  a  very  short  distance 
in  ice.  We  were  able  to  demonstrate  conclusively  that  the 
laser  beam  can  cut  cleanly  and  rapidly  through  both  firn  and 
ice  and  that  it  can  be  manipulated  efficiently  with  standard 
optical  systems. 


Subsequently,  upon  our  return  to  the  laboratory  at  the 
University  of  Kansas,  we  tested  the  system  on  a  4  inch  diameter 
ice  core  from  the  Greenland  ice  sheet  taken  at  a  depth  of  about 
170  meters  (Koci,  personal  communication).  The  core  was 
supplied  to  us  by  P.  Grootes.  in  this  case,  cutting  was 
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performed  in  an  open  freezer  with  the  initial  temperature  of 
approximately  -20  C  using  about  15  watts  total  power  output. 

The  cutting  was  performed  without  optical  condensers  using  the 
beam  directly  from  the  laser.  We  were  able  to  slice  off  a  2 
inch  segment  of  the  core  without  shattering  or  fragmentation  of 
any  kind.  The  cut  was  about  2  mm  wide  and  could  have  been 
reduced  to  half  that  width  with  an  optical  system.  A  video 
tape  taken  in  the  laboratory  at  the  time  of  the  cutting 
experiment  was  prepared  and  sent  to  NSF/DPP  and  PICO. 

In  the  process  of  completing  this  experiment,  we 
determined  that  it  would  be  possible  to  develop  an  optical 
system  that  would  permit  the  beam  to  be  rotated  in  a  circular 
path  that  could  be  used  for  cutting  deep  ice  cores.  With  minor 
modifications  this  system  could  be  used  in  fluid  filled  holes 
as  well  as  in  open  holes. 

The  fact  that  the  laser  beam  cuts  entirely  by  melting  and 
exerts  no  torque  on  the  ice  at  the  cutting  surface,  greatly 
reduces  the  potential  for  fracturing  and  breakage  of  the  core 
during  drilling.  Even  more  important,  the  fact  that  the  beam 
can  be  deflected  at  90  degrees  to  the  drilling  direction  means 
that  the  ice  core  can  be  cut  off  at  the  bottom  of  the  hole  by 
the  beam.  This  has  the  effect  of  greatly  reducing  stress  on 
the  core  at  the  time  it  is  lifted  free  of  the  bottom  of  the 
hole.  Finally,  chemical  contamination  cannot  be  introduced  by 
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the  laser  beam  and  the  infrared  wavelengths  used  for  drilling 
are  too  long  to  cause  significant  radiolytic  breakdown  keeping 
chemical  alterations  of  all  types  at  an  absolute  minimum. 

In  construction  of  a  prototype,  the  rotating  mass  of  the 
optical  system  and  the  thin-walled  core  barrel  will  be  so  low 
that  only  the  lightest  of  anti -torque  systems  must  be  used. 

The  entire  drill  assembly  including  the  laser,  optics,  core 
barrel  and  scavenger  pump  system  can  be  expected  to  weigh  less 
than  200  lbs  and  can  be  supported  on  a  light  weight  cable  that 
must  include  electrical  conductors  for  the  28  volt  power  supply 
to  the  laser.  Power  requirements  are  modest.  Electrical  power 
to  the  completed  system  can  be  supplied  by  a  3,000  watt 
generating  facility  and  should  be  fully  adequate  for  all 
operations  Including  drilling  and  hoisting  of  the  drill  and 
core  barrel  assembly. 

The  simplest  drill  design  is  shown  in  Figure  A.  In  this 
design,  the  carbon  dioxide  laser  is  mounted  vertically  in  the 
bore  hole  and  a  rotating  head  containing  deflecting  mirrors  is 
attached  so  that  the  laser  beam  can  be  directed  down  against 
the  ice  and  turned  in  a  circle  by  a  small,  low  power,  direct 
current  motor.  In  principle,  it  is  possible  to  design  the 
drill  so  that  the  radius  of  the  core  is  adjustable  but  we 
propose  to  build  a  system  with  a  fixed  core  diameter  (most 
probably  4  inches ) . 
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In  any  vertical  ice  drilling  operation  involving  an 
infrared  laser  beam  it  is  necessary  to  remove  the  water  from 
the  annulus  cut  by  the  beam.  This  will  be  done  by  a  small 
scavenger  piunp  that  will  pick  up  the  water  at  the  point  where 
it  is  being  produced  and  transfer  it  to  a  holding  tank  at  the 
top  of  the  core  barrel.  It  will  be  necessary  to  empty  this 
holding  tank  each  time  the  core  barrel  is  withdrawn  from  the 
hole. 

Figure  B,  is  a  design  for  an  optical  system  that  would 
permit  switching  the  laser  beam  path  from  the  position  used  for 
drilling  (through  beam  window  1)  to  the  position  used  for 
cutting  the  core  lose  from  the  bottom  of  the  hole  (through  beeun 
window  2).  This  would  be  accomplished  by  activating  the  beam 
selector  at  the  center  of  the  optical  bridge.  This  portion  of 
the  drill  assembly  would  be  located  at  the  top  of  the  core 
barrel . 

Figure  C,  shows  the  design  of  the  core  barrel.  It  would  be 
constructed  of  two  thin-wall  steel  tubes  arranged 
concentrically  with  the  space  between  the  tubes  being  used  as 
the  location  of  the  optical  wave  guide  fiber  and  the  vacuum 
line  to  the  scavenger  pump.  In  the  proposed  prototype  drill, 
we  will  construct  a  core  barrel  to  accomodate  a  1.5  meter  ice 
core  segment  to  reduce  weight  and  make  preliminary  testing 
easier.  In  a  scaled  up  version,  the  core  barrel  could  be 
extended  to  3  or  4  meters. 
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Using  the  laser  that  we  already  have  and  some  of  the 
equipment  from  our  laboratory,  we  expect  to  produce  a  proof-of- 
concept  prototype.  Because  of  the  ruggedness  and  high  efficiency 
of  carbon  dioxide  lasers,  we  believe  that  this  system  could  be 
upgraded  to  produce  a  workable  drill  for  use  in  polar  regions 
that  would  greatly  reduce  the  logistic  requirements  for  transport 
and  power. 
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FIGURE  CAPTIONS 


Figure  A.  The  drill  assembly  consists  mainly  o-f  an  infra-red 
C02  laser,  optics,  core  barrel  and  scavenger  pump  system,  and 
is  e::pected  to  weigh  less  than  200  lbs. 

Figure  B.  Design  of  the  optical  system  that  would  permit 
switching  the  laser  beam  path  from  the  position  for  drilling 
to  the  position  for  cutting  the  core  loose  at  the  bottom  of 
the  drill  hole. 

Figure  C.  The  design  of  the  core  barrel  includes  the 
construction  of  two  thir*— walled  concentric  steel  tubes.  The 
space  between  the  tubes  will  be  used  as  tlie  location  of  the 
wave  guide  fiber  and  the  vacuum  line  to  the  scavenger'  pump. 
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(as  they  predicted,  given  the  complex  profiles);  this  can  be 
explained  by  a  precessional  variation  of  unresolved  components 
of  the  He  ii  profile  in  their  analysis. 

The  corresponding  masses  of  the  compact  object  and  the 
normal  star  can  be  easily  derived  from  equation  1;  we  find 
M*  =  3.2±0.4Mo  and  M,==0.8±0.1M®.  The  above  estimates 
for  F  and  q,  taken  at  face  value,  imply  therefore  that  the  compact 
object  is  a  neutron  star  rather  than  a  black  hole.  Given  the  high 
quality  of  our  spectra,  we  are  fairly  confident  of  the  value  of 
K.  A  somewhat  higher  mass  ratio  could  be  inferred  considering, 
as  suggested  in  ref.  16,  a  moderate  flattening  of  the  secondary 
or  a  wind-type  mass  flow  that  would  obscure  the  X-rays  for  a 
small  fraction  of  the  eclipse. 

The  quality  of  our  data  also  allows  us  to  study  the  position 
and  strength  of  the  line  peaks.  For  this  purpose  we  chose  seven 
spectra  with  clear  double-peaked  structure.  Note  that  the  peaks 
are  rather  weak  (Fig.  1)  so  that  the  line  always  shows  a  strong 
central  component  with  the  orbital  variation  described  above. 
Sine  fits  to  the  position  of  the  two  underlying  gaussians  and  to 
the  intensity  ratio  of  the  two  peaks  are  given  in  Table  1.  We 
find  (Fig.  2)  that  the  red  peak  is  stronger  when  the  line  is 
redshifted  and  the  blue  peak  is  stronger  when  the  line  is  blue- 
shifted.  The  position  of  the  blue  peak  seems  to  vary  much  less 
than  the  position  of  the  red  peak. 

Unequal  orbital  modulation  of  the  two  peaks  seems  to  support 
the  idea  of  two  unrelated  emission  regions,  one  approaching  us 
and  the  other  receding  at  several  hundred  kilometres  per  second. 
But  the  He  II  line  at  4,686  A  is  produced  in  a  hot  environment, 
so  one  of  the  peaks  is  likely  to  be  emitted  in  the  binary  system 
region  which  should  have  small  systemic  velocity. 

In  cataclysmic  vari.ibles  as  well  as  X-ray  binaries,  the  He  ii 
emission  comes  from  the  vicinity  of  the  compact  object  or,  more 
precisely,  from  the  corona  surrounding  the  accretion  disk.  We 
believe  that  this  is  the  case  in  SS433  as  well.  The  peak  separation 
of  ~7.8  A  corresponds  to  the  keplerian  radius  of  «  3  x  10"  cm 
for  a  neutron  star  that  is  smaller  than  the  radius  of  its  Roche  lobe. 

The  intensity  ratio  of  the  two  peaks  is  such  that  the  part  of 
the  accretion  disk  facing  the  normal  star  is  always  dimmer.  This 
rules  out  the  illumination  from  the  normal  star  as  a  possible 
source  of  different  peak  strengths.  Similarly,  the  eclipse  of  the 
accretion  disk  is  far  too  short  to  cause  this  effect.  Another 
possibility  would  be  the  contribution  by  emission  from  the  hot 
spot;  it  is  unlikely,  however,  that  the  hot  spot  lies  on  the  part 
of  the  accretion  disk  facing  away  from  the  secondary.  Finally, 
the  emission  of  the  He  ii  line  from  the  part  of  the  corona  facing 
the  donor  star  might  be  partly  absorbed  by  a  strong  wind  from 
the  hot-spot  region.  The  existence  of  such  outflow  has  been 
inferred  from  the  analysis  of  photometric  light  curves"  and 
X-ray  observations'*.  The  rotation  of  the  disc  together  with  such 
periodic  obscuration  can  explain  both  the  strength  and  the 
separation  modulation  of  the  peaks.  On  the  other  hand,  the  high 
scatter  of  the  redshifted  line  modulation  and  the  large  separation 
of  the  peaks  on  the  spectrum  obtained  on  18  October  indicate 
additional  nonperiodic  effects.  Double  peaks  have  also  been 
observed  in  the  Paschen  stationary  lines'’.  The  splitting  is  smal¬ 
ler  than  observed  in  the  He  it  line.  This  is  consistent  with  the 
picture  of  lower-excitation  lines  being  emitted  in  the  outer  parts 
of  the  disk  (which  are  cooler  and  rotate  more  slowly)  or  in  the 
wind  from  it. 

In  conclusion,  the  lower  value  of  K  measured  here,  together 
with  the  estimate  of  the  mass  ratio  from  X-ray  observations  of 
the  primary  eclipse,  lead  us  to  believe  that  the  compact  object 
in  SS433  has  a  lower  mass  than  previously  thought.  The  double¬ 
peak  structure  observed  in  the  He  ti  line  at  4,686  A  is  here 
interpreted  as  due  to  emission  from  the  corona  of  the  accretion 
disk.  □ 
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Detection  of  the  small  flux  of  extraterrestrial  neutrinos  expec¬ 
ted  at  energies  above  1  TeV,  and  identification  of  their  astro- 
physicai  point  sources,  will  require  neutrino  telescopes  with 
effective  areas  measured  in  square  kilometres — much  larger  than 
detectors  now  existing''^.  Such  a  device  can  be  built  only  by  using 
some  naturally  occurring  detecting  medium  of  enormous  extent: 
deep  Antarctic  ice  is  a  strong  candidate.  A  neutrino  telescope  could 
be  constructed  by  drilling  holes  in  the  ice  with  hot  water  into  which 
photomultiplier  tubes  could  be  placed  to  a  depth  of  1  km.  Neutrinos 
would  be  recorded,  as  in  underground  neutrino  detectors  using 
water  as  the  medium,  by  the  observation  of  Cerenkov  radiation 
from  secondary  muons.  We  have  begun  the  AMANDA  (Antarctic 
Muon  and  Neutrino  Detector  Array)  project  to  test  this  idea,  and 
here  we  describe  a  pilot  experiment  using  photomultiplier  tubes 
placed  into  Arctic  ice  in  Greenland.  Cerenkov  radiation  from 
muons  was  detected,  and  a  comparison  of  count  rate  with  the 
expected  muon  flux  indicates  that  the  ice  is  very  transparent,  with 
an  absorption  length  greater  than  18  m.  Our  results  suggest  that 
a  full-scale  Antarctic  ice  detector  is  technically  quite  feasible. 

Neutrinos  would  be  detected  by  the  standard  technique  of 
observing  the  Cerenkov  light  emitted  by  upward-moving  muons, 
produced  by  neutrinos  interacting  in  or  below  the  detector.  At 
these  high  energies,  the  direction  of  the  muon  produced  is  within 
1°  of  the  direction  of  the  parent  neutrino,  so  the  observation  of 
astrophysical  neutrino  sources  is  possible  in  principle.  Such  a 
detector  at  the  South  Pole  would  have  many  advantages:  (1) 
unlike  existing  and  some  proposed  detectors,  there  would  be 
no  physical  limit  to  its  size  and  structure,  and  it  could  be 
gradually  expanded  (holes  could  be  added  year  after  year);  (2) 
the  ice  would  form  a  mechanical  structure  for  the  instrument, 
enabling  all  active  electronics  other  than  photomultiplier  tubes 
(PMTs)  to  be  placed  directly  above  at  the  surface,  putting  most 
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of  the  electronics  within  reach  for  adjustment  and  maintenance; 
(3)  the  lack  of  background  light  sources  (for  example,  natural 
radioactivity)  deep  in  the  ice  combined  with  low  ambient  tem¬ 
peratures  imply  that  PMT  noise  rates  should  be  very  low;  (4) 
location  at  the  South  Pole  would  have  unique  astronomical 
advantages — for  instance,  the  observation  of  sources  at  constant 
zenith  angle  and  the  ability  to  observe  neutrinos  from  a  Northern 
Hemisphere  source  simultaneously  with  the  observation  of  ultra- 
high-energy  gamma  rays  from  the  same  source  by  Northern 
Hemisphere  air-shower  arrays.  For  this  technique  to  work,  the 
ice  must  be  very  transparent  (peak  optical  absorption  length  > 
several  metres).  In  this  initial  study,  conducted  at  the  Greenland 
Ice  Sheet  Project  II  site  (72°  N,  38°  W.)  during  13-16  August 
1990,  muon  rate  measurements  were  made  to  assess  the  optical 
clarity  of  polar  ice. 

A  string  of  three  Thom  EMI  model  9870  PMTs,  each  12  cm 
in  diameter  with  a  hemispherical  window,  was  lowered  down  a 
IS-cm-diameter  borehole  to  a  depth  of  -217  m,  ~100  m  below 
the  hm  (packed  snow)  layer  (Fig.  1).  Operating  close  to  single 
photoelectron  sensitivity,  the  PMTs  had  noise  rates  between  10 
and  20  kHz,  higher  than  expected,  and  due  to  surface  light 
leaking  down  the  borehole.  The  event  trigger  required  that  all 
three  PMTs  produce  a  digital  pulse  within  30  ns  of  each  other; 
the  observed  coincidence  rate  was  1.8  Hz.  By  inserting  additional 
time  delays  to  destroy  the  real  coincidence  window,  the  purely 
random  coincidence  rate  was  directly  measured  to  be  0.2  Hz, 
well  below  the  1.8-Hz  coincidence  rate.  From  previous  muon 
flux  measurements  and  calculations'*,  the  vertical  muon  flux  at 
depth  (217  m  =  180  m  water-equivalent)  is  1  muon  m'^  s”‘  sr"*. 
The  corrected  event  rate  of  1.6  Hz,  together  with  the  expected 
vertical  flux,  indicates  that  each  PMT  was  sensitive  to  muons 
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FIG.  1  PMT  string  used  in  Greenland.  The  three  PMT s  were  arranged  vertically, 
with  1-m  separation  between  PMTs.  Each  PMT  had  a  hemispherical  window, 
and  faced  upwards.  The  hole  was  filled  with  butyl  acetate,  an  organic  liquid 
chosen  for  its  low  freezing  point  and  clarity,  to  a  depth  of  11  m.  for  good 
optical  contact  between  the  PMTs  and  the  borehole  walls.  The  PMTs  were 
equipped  with  voltage  divider  chains  housed  in  pressure  vessels,  and  high 
voltage  was  provided  from  the  surface  directly.  The  PMT  analog  output  was 
routed  through  a  coaxial  cable  (RG  180)  directly  to  a  receiver  on  the  surface 
that  was  specifically  designed  by  Phillips  Electronics  to  compensate  for  the 
attenuation  and  dispersion  introduced  by  the  cable.  Using  this  receiver,  only 
a  small  attenuation  (<25%)  was  observed  in  the  amplitude  of  the  PMT 
output,  and  the  characteristic  ~10ns  rise  time  of  the  PMT  signal  was 
preserved.  The  output  signals  from  the  receiver  were  routed  into  amplifiers, 
coincidence  circuitry  and  CAMAC  readout  electronics.  Discriminator  thresh¬ 
olds  were  set  as  low  as  possible  (30  mV),  resulting  in  a  sensitivity  to  signals 
greater  than  1  photoelectron  equivalent  (p.e.)  for  the  top  PMT,  1  p.e.  for  the 
middle  PMT  and  2  p.e.  for  the  bottom  P^. 


arriving  within  a  radius  of  ~1  m  (Fig.  2).  Thus,  the  experiment 
shows  that  a  PMT  placed  at  a  depth  of  200  m  in  polar  ice  is 
capable  of  acting  as  a  detector  with  an  effective  area  >3  m^.  In 
addition,  the  Monte  Carlo  studies  described  next  indicate  that 
if  all  PMTs  were  operated  at  the  single  photoelectron  level  (not 
possible  in  Greenland  because  of  the  high  noise  rates),  the 
effective  area  of  the  string  would  increase  from  ~3  to  —6  m^. 

Ice  at  this  depth  has  a  large  concentration  of  air  bubbles;  the 
ice  cores  taken  from  this  depth  are  milky  and  translucent.  Davis’ 
has  calculated  the  non-isotropic  intensity  distribution  of  light 
scattered  by  air  bubbles  in  water,  which  is  a  good  approximation 
to  the  case  of  air  bubbles  in  ice.  We  have  done  Monte  Carlo 
calculations  of  photon  propagation  in  the  ice,  using  the  Davis 
results.  In  the  limit  where  Ns 20,  a  numerical  fit  to  the  radial 
distributions  generated  by  the  Monte  Carlo  gives  d  = 
A(1  -I•3■^/N),  where  d  is  the  mean  radial  distance  reached  after 
N  scatters,  and  A  is  the  mean  free  path  for  light  scattering  and 
A  =  l/(crp),  where  cr  is  the  bubble  cross-sectional  area  and  p  is 
the  number  density  of  bubbles.  Analysis  of  cores  taken  pre¬ 
viously  from  this  depth  in  Greenland  (A.  J.  Gow,  personal 
communication)  indicates  that  there  are  200  bubbles  cm"’  of 
mean  radius  0.017  cm,  giving  A  ~ 5.5  cm.  (For  our  particular 
experiment,  N  =  40,  using  d  =  1  m  from  the  rate  analysis  given 
previously  and  A  =  5.5  cm).  These  results  have  been  incorporated 
into  a  Monte  Carlo  calculation  to  predict  muon  rates  and  pulse- 
height  distributions  for  each  PMT.  Using  the  known  muon  flux 
and  angular  distribution,  Cerenkov  photons  were  generated 
from  the  muon  paths  and  propagated  through  the  bubbly  ice. 
The  program  included  the  quantum  efficiency  and  collection 
efficiency  of  the  PMTs  as  measured  by  the  manufacturer  and 
the  variation  of  PMT  effective  area  with  the  zenith  angle  of  an 
incoming  photon,  measured  as  part  of  the  calibration  of  this 
type  of  PMT  for  the  Irvine-Michigan-Brookhaven  proton  decay 
experiment’.  The  absorption  length  of  the  light  in  the  clear  ice 
between  bubbles  was  scaled  as  a  fraction  of  measured  laboratory 
values’  (it  should  be  pointed  out  that  laboratory  measurements 
of  absorption  include  some  scattering  effects,  and  may  therefore 
underestimate  the  true  absorption  length).  Good  agreement 
between  the  predicted  and  observed  pulse-height  distributions 
is  achieved  only  if  the  peak  optical  absorption  length  of  the  in 
situ  polar  ice  is  greater  than  ~18  m,  or  75%  of  the  absorption 
length  measured  for  pure  laboratory  ice  (Fig.  3).  Systematic 
uncertainties  in  the  Monte  Carlo  calculation  (for  example,  an 


FIG.  2  The  PMT  string  shown  as  a  cylindrical  detector,  for  calculation  from 
the  muon  rate  of  the  effective  area  of  a  PMT.  The  geometrical  aperture  of 
the  instrument,  in  m^  sr.  times  the  muon  flux  (1.0  m'^  sr"‘  s"‘)  yields  the 
muon  rate  (1.6  Hz),  so  the  aperture  is  1,6  m^  sr.  Visualizing  the  F^  string 
as  shown,  the  geometrical  aperture  is** 


-j  (2Rf  „  +  h*  -  V4f?^„h’+h*) 


This  yields  R,„=1.0m,  or  A,„=3.1  m^. 
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FIG.  3  Integral  muon  rate  as  a  function  of  minimum  pulse  height  (in  p.e.) 
required  of  all  three  PMTs.  Filled  circles  show  experimental  results,  solid 
line  shows  Monte  Carlo  predictions  If  the  Greenland  ice  has  peak  optical 
absorption  length  of  25  m  (identical  to  laboratory  ice).  Dotted  lines  show 
Monte  Carlo  predictions  for  peak  absorption  lengths  of  18  m  and  12  m.  The 
error  bars  show  only  statistical  uncertainties  in  the  observed  rates  and  in 
the  Monte  Carlo  rates. 


observed  presence  of  impurities  in  the  butyl  acetate  used)  would 
tend  to  reduce  the  predicted  Monte  Carlo  rates,  so  our  results 
should  be  regarded  as  a  lower  limit  to  the  optical  clarity  of  the 
ice.  It  should  also  be  pointed  out  that  the  large  absorption  length 
implied  by  these  results  indicates  that  PMT  size,  rather  than  the 
optical  quality  of  the  ice,  will  be  the  limiting  factor  for  a  real 
neutrino  detector;  this  is  because  for  a  smaller  PMT,  loss  of 
light  due  to  angular  spread  of  the  Cerenkov  cone  outweighs  the 
loss  of  light  due  to  optical  absorpf  ion  of  he  ice.  Finally,  because 
the  ionic  purity  and  particulate  c<.  ncen  .ation  of  the  Greenland 
ice  are  comparable  to  those  measured  for  distilled  water,  it  is 
not  surprising  that  the  optical  properties  of  the  two  substances 
are  also  comparable.  Trodahl  and  others’  have  studied  the 
diffusive  transport  of  light  through  bubbly  sea  ice  and  conclude 
that  the  optical  absorption  length  at  a  wavelength  of  SOO  nm 
must  be  at  least  10  m;  as  this  ice  contains  brine  pockets  and 
algae,  the  pure  ice  found  in  the  interior  of  Greenland  or  Antarc¬ 
tica  should  be  much  better  for  our  purposes. 

We  find  these  results  very  encouraging,  and  are  planning  more 
extensive  experiments  at  the  South  Pole  during  the  coming 
austral  summer  (November  1991-January  1992),  in  particular 
the  operation  of  one  or  two  strings,  each  with  four  20-cm- 
diameter  PMTs  at  a  depth  of  ~  1  km,  where  the  ice  is  expected 
to  be  bubble  free^  Assuming  that  these  tests  are  successful,  we 
then  hope  to  begin  construction  and  installation  of  full-size 
strings  of  twenty  20-cm-diameter  PMTs  in  the  succeeding  season 
(1992-93).  In  addition,  as  part  of  this  year’s  work,  we  will  be 
investigating  the  possibility  of  using  our  detection  technique  to 
construct  muon  detectors  able  to  operate  in  conjunction  with 
South  Pole  Air  Shower  Experiment,  an  existing  air-shower  array 
at  the  South  Pole'®.  The  large  effective  muon-sensitive  area  of 
a  PMT  at  shallower  depths  suggests  that  good  muon  coverage 
could  be  achieved  with  only  a  small  number  of  drilled  holes.  □ 
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The  discovery  of  superconductivity  in  potassium-doped  C«o(ref.  1) 
has  been  followed  by  an  intense  effort  to  understand  the  physics 
and  chemistry  of  metal-doped  fullerene  solids’'”.  Experimental 
studies  of  alkali-metal-doped  C«o  have  now  provided  insight  into 
the  structure^*'’  and  the  coherence  length  and  penetration  depth’ 
of  the  superconducting  phase.  No  measurements  of  the  super¬ 
conducting  energy  gap  (A)  have,  however,  been  reported.  The  BCS 
theory  of  superconductivity”,  which  has  been  used  to  interpret 
much  of  this  experimental  work’’’’*''’,  predicts  (in  the  limit  of 
weak  coupling)  that  the  reduced  energy  gap  has  a  material- 

independent  value  of  3.53.  Values  in  excess  of  3.5  define  strong 
coupling,  and  thus  provide  insight  into  the  nature  of  the  pairing 
mechanism.  Here  we  describe  the  measurement  of  A  for  single¬ 
phase  superconducting  Rb}C«o  by  tunnelling  spectroscopy  using  a 
scanning  tunnelling  microscope.  We  obtain  a  value  of  A  at  4.2  K 
of  6.6  ±0.4  meV,  corresponding  to  a  reduced  energy  gap  of  53. 
This  is  significantly  larger  than  predicted  by  BCS  theory,  but 
similar  in  magnitude  to  values  found  for  high-temperature  copper 
oxide  superconductors”.  Our  finding  of  strong  coupling  in  Rb3CM 
suggests  the  need  for  caution  in  using  standard  BCS  theory  to 
interpret  superconductivity  in  metal-doped  Cto- 

Tunnelling  spectroscopy  has  been  one  of  the  most  successful 
techniques  used  to  probe  conventional  metal  and  alloy  super¬ 
conductors'*.  In  particular,  tunnelling  spectra  can  provide 
values  for  both  the  energy  gap  A  and  the  electron-phonon 
spectral  function  a’F(<u),  where  a’  is  a  measure  of  the  coupling 
and  F(ii>)  is  the  phonon  density  of  states.  To  obtain  a  clear 
value  for  A  by  tunnelling,  it  is  essential  to  prepare  a  uniform 
insulating  barrier  between  the  superconductor  and  metal  junc¬ 
tion.  As  the  Cjo  superconductors  have  short  coherence  lengths’ 
and  as  the  samples  are  inhomogenous,  we  expect  that  planar 
junctions  might  show  significantly  broadened  gap  features  (even 
for  the  ideal  BCS  case).  Wc  have  therefore  used  a  low- 
temperature  scanning  tunnelling  microscope  (STM)  to  make 
point  junctions  with  a  sharpened  metal  tip.  Point  tunnelling 
spectroscopy  has  been  used  to  measure  A  in  a  variety  of  super¬ 
conducting  materials,  including  pure  metals'*,  organics'’  and 
oxides”,  and  the  values  of  A  determined  by  this  technique  for 
conventional  superconductors  are  comparable  to  those  obtained 
with  planar  junctions. 

The  Rb3C#o  samples  were  prepared  by  methods  described  in 
detail  elsewhere'’.  Briefly,  purified  and  dried  Cm  was  reacted 
with  RbHg  in  a  1 : 3  ratio  at  200  °C  for  8-72  h.  After  this  initial 
reaction,  the  superconducting  fraction  determined  from  shield¬ 
ing  measurements  was  typically  >35%.  The  resulting  powder 
was  then  ground,  pressed  into  a  pellet,  and  sintered  at  200 ‘C 
for  3-12  h.  Diamagnetic  shielding  measurements  made  on  these 
sintered  pellets  show  that  the  superconducting  fraction  can 


*  To  wtwfn  oorrotDOodtnco  sPoiAd  bo  odorooooA. 
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